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ABSTRACT 

We present results from three-dimensional hydrodynamic simulations of the high redshift collapse of 
pregalactic clouds including feedback effects from a soft H2 photodissociating UV radiation field. The 
simulations use an Eulerian adaptive mesh refinement technique to follow the nonequilibrium chemistry 
of nine chemical species with cosmological initial conditions drawn from a popular A-dominated cold 
dark matter model. The results confirm that the soft UV background can delay the cooling and collapse 
of small halos (~ IO^Mq). For reasonable values of the photo-dissociating flux, the H2 fraction is in 
equilibrium throughout most of the objects we simulate. We determine the mass threshold for collapse 
for a range of soft-UV fluxes and also derive a simple analytic expression. Continuing the simulations 
beyond the point of initial collapse demonstrates that the fraction of gas which can cool depends mostly 
on the virial mass of the halo and the amount of soft-UV flux, with remarkably little scatter. We 
parameterize this relation, for use in semi-analytic models. 

Subject headings: cosmology: theory - early universe - galaxies: formation 



1. INTRODUCTION 

A critical epoch in the evolution of the universe is the 
the post-recombination period during which the first stars 
form, affect their environment and help determine how the 
intergalactic medium (IGM) ultimately becomes reionized. 
Recent observations are beginning to constrain models of 
this epoch. The observation of quasars at redshifts as high 
as z = 5.8 (Fan et al. 2000) with no accompanying Gunn- 
Peterson absorption trough in their spectra indicates that 
the onset of reionization occurs quite early {z > 7) and 
so the first generation of stars forms even earlier. Ob- 
servations of quasar spectra (Ellison et al. 1999, Ellison 
et al. 2000) suggest that metallicities possibly as high as 
logCIV/ HI — —2.6 may be present even in the low col- 
umn density log A^hi < 14.5 Lyman alpha forest. This too 
may be indicative of an active period of pre enrichment 
of the IGM by an early generation of stars (Ostrikcr & 
Gnedin 1996). Future multiwavelength observations may 
be sensitive probes of this star formation epoch and the 
subsequent onset of reionization. A partial list of pro- 
posed observations include direct imaging of quasars and 
star clusters at z > 10 by the Next Generation Space Tele- 
scope (Haiman & Loeb 1999, Barkana & Loeb 2000), ab- 
sorption measurements in the 21 cm line using the Square 
Kilometer Array to identify the first epoch of massive star 
formation (Tozzi et al. 1999), and searches for secondary 
anisotropics in the Cosmic Microwave Background by the 
next generation of millimeter telescopes such as ALMA 
to constrain properties of the reionization epoch (Benson 
et al. 2000). Thus it is important to understand what our 
best cosmological models predict for the evolution of these 
early star forming structures and what governs the amount 



of gas that can ultimately condense into stars. 

It has long been understood that a high redshift popula- 
tion of stars could form naturally in hierarchical models of 
structure formation. In these models small masses collapse 
first and later merge into larger structures. The ability of 
these objects to form stars depends on the ability of gas 
within these objects to cool radiatively, reduce pressure 
support, and condense. However, the first structures to 
become self-gravitating in these models would have tem- 
peratures T < 10^ K , too cool for HI line cooling to be 
effective. The only other coolant available to gas with a 
primordial composition is molecular hydrogen. Therefore 
it is the complex chemistry of H2 in these objects that 
controls whether the gas can condense to sufficiently high 
densities for star formation to take place. 

Early analytical models of the collapse of the first small 
scale structures studied the effect of H2 cooling within 
highly idealized cloud geometries, typically a homogeneous 
gas cloud radiatively cooling while undergoing a uniform 
free fall collapse (Saslow & Zipoy 1967, Peebles & Dicke 
1968, Hirasawa 1969, Matsudo et al. 1969, Hutchins 1976, 
Silk 1983, Palla et al. 1983, Lepp & ShuU 1984, Lahov 
1986). 1-D simulations sought to improve upon these an- 
alytical estimates, but still had to impose the shape of 
the initial density profiles and mass of the virialized ob- 
ject as initial conditions within an assumed cloud geome- 
try (Haiman, Thoul & Loeb 1996, Tegmark et al. 1997, 
Omukai & Nishi 1998, Nishi et al. 1998, Nakamura & 
Umemura 1999). Haiman, Thoul & Loeb (1996) followed 
the dynamics of gas and cold dark matter and the H2 
chemistry in a single, spherical density peak and showed 
the importance of shell crossing between the gas and dark 



1 Hubble Fellow 



2 



matter in determining the fraction of mass within an ob- 
ject that could coUapse and the minimum mass for col- 
lapse. In the same spirit Tegmark et al. (1997) used top 
hat density profiles and a simplified collapse criterion, i.e. 
that the cooling time be less than the Hubble time for the 
object right after virialization, to obtain rough estimates 
of the critical H2 mass fraction needed for cooling and the 
minimum mass (as a function of virialization redshift) that 
could produce it. However, the minimum mass to collapse 
determined in this way is strongly dependent on the H2 
cooling function used. Abel et al. (1998) and Fuller & 
Couchman (2000) recalculated this minimum mass for top 
hat density profiles using the Tegmark et al. (1997) col- 
lapse criterion with updated H2 cooling functions (Lepp & 
Shull 1983 and Galh & Palla 1998, respectively) and found 
that the minimum mass able to collapse (at virialization 
rcdshifts ~ 30) was reduced by as much as an order of 
magnitude. 

Molecular hydrogen is fragile and easily photodissoci- 
ated by photons whose energies lie below 13.6 eV (Stecker 
& Williams 1967). Thus even before reionization, radia- 
tion from the first stars could affect H2 abundances at 
large distances (Dekel & Rees 1987). Haiman, Rees, & 
Loeb (1997) quantified this negative radiative feedback in 
more detail. They studied the collapse of homogeneous, 
spherical gas clouds in the presence of a soft UV photodis- 
sociating flux assuming both ionization and chemical equi- 
librium within the cloud to determine which masses could 
still collapse. Haiman, Abel & Rees (1999) improved on 
that calculation by solving spherically symmetric radiative 
transfer coupled with the time dependent rate equations 
for the H2 chemistry across truncated isothermal spheres. 
They also modified the collapse criterion such that a struc- 
ture was said to be able to collapse and form stars if at 
any time during the object's history the H2 cooling time in 
the core became less than the time elapsed since the core's 
formation. This allowed the fraction of H2 abundance re- 
quired for cooling and collapse of the core to be lower in 
certain cases than previous estimates of the "critical" H2 
fraction based on the Tegmark et al. (1997) collapse crite- 
rion. Analytical models that include a first population of 
OB-like stars have been used to probe the suppression of 
star formation elsewhere within the host cloud due to the 
destruction of H2 coolant by the direct radiation from the 
first stars (Omukai & Nishi 1999; Glover & Brand 2000). 

The major drawback affecting all of the above studies 
whether with or without radiative feedback is their one 
dimensional nature. Hierarchical structure formation that 
predicts the growth of structure in and along filaments is 
inherently three dimensional. However, fully three dimen- 
sional simulations are computationally challenging due to 
the need for large dynamic range to model simultaneously 
the internal structure of the first star forming structures 
and the effects of gravitational tidal forces and merging on 
their growth. Such 3-D simulations have only recently be- 
come computationally feasible. Gnedin & Ostriker (1997) 
simulated the thermal history of the intergalactic medium 
and the population of first objects to collapse but had in- 
sufficient resolution to probe inside these first structures. 
Abel et al. (1998) used a static nested grid approach 
and constrained realizations of 3 ~ 4(7 primordial fluc- 
tuations in a cluster normalized SCDM model to study 
when these structures were able to cool, but were also lim- 



ited by spatial resolution 1 kpc) from probing deeply 
within the peaks. Subsequent higher resolution 3-D sim- 
ulations focused on the properties of a single collapsing 
cloud. Bromm, Coppi, & Larson (1999) studied the frag- 
mentation of a 2 X 1O^M0 rotating disk under highly ide- 
alized initial conditions using a TREESPH code. Fuller 
& Couchman (2000) used an adaptive smoothing-kernal 
SPH coupled to a P'^M gravity solver in a simulation box 
of comoving length 25h'~^ kpc to investigate the proper- 
ties of the first peak to collapse. However, because of 
the small box size they needed to impose an overdense 
cold dark matter cosmology to model the formation of 
the cloud within its host filament. The highest spatial 
and mass resolution self-consistent cosmological 3-D sim- 
ulations to date of the collapse and fragmentation of pri- 
mordial clouds (Abel, Bryan, & Norman 1998, 2000; Nor- 
man, Abel, & Bryan 2000) use Eulerian adaptive mesh 
refinement (AMR) techniques to achieve extremely large 
dynamic range (currently > 3 x 10^). They start from cos- 
mological initial conditions in a cluster normalized SCDM 
model to study the detailed internal evolution of a single, 
collapsing cloud in interaction with its environment in a 
box with comoving length L = 128 kpc. 

In this paper we present the results of fully 3-D Eu- 
lerian AMR simulations of the formation and collapse of 
primordial pregalactic structure in a realistic ACDM cos- 
mology. This study expands upon previous work in several 
important ways. First, we use a larger simulation volume 
(1 Mpc^ comoving) to model more accurately the gravita- 
tional tidal forces and effects of merging on the evolving 
structures. Second, we introduce a simple model for ra- 
diative feedback. We analyze the collapse properties of in- 
dividual peaks with and without the presence of this soft 
UV background and compare our results with previous 
1-D studies. Finally, we use the population of pregalac- 
tic objects evolving together in the most active structure 
forming region of our simulation to develop statistics on 
the fraction of gas that is able to cool due to the presence 
of H2 and the fraction of gas that is both cold and dense 
enough to be available for star formation as functions of 
cloud mass and the average soft UV radiation field it ex- 
periences. These statistics represent the key new result 
of this paper. The amount of gas available for star for- 
mation in pregalactic objects is a critical input parameter 
into recent semi-analytical attempts to model stellar feed- 
back from the first stars on the intergalactic medium, sub- 
sequent structure formation, and the nature of reioniza- 
tion (Gnedin & Ostriker 1997, Ciardi et al. 2000, Ferrara 
et al. 2000, Ciardi, Ferrara, & Abel 2000). 

While we focus on the negative effect of radiation in this 
paper, it should be noted that positive feedback is possi- 
ble if the first objects also emit high energy X-ray photons 
(Haiman, Abel & Rees 1999; Oh 2000). We wifl address 
this point in future work. 

This paper is organized in the following way. In Sec- 
tion §^ we describe the set-up of our simulations including 
simulation parameters and technique, cosmological model, 
H2 chemistry, and radiative feedback. In Section §^ we 
present the characteristics of the simulated data set of pre- 
galactic objects identified in our simulations and used for 
further analysis. In Section ^ we use the data set as a 
whole to develop fitting formulaes for the fraction of gas 
that can cool, become dense and thus be available for star 
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formation with and without the presence of a background 
H2 photodissociating flux. In Section §|| we use radial pro- 
files to show the internal collapse characteristics of these 
structures. We consider both the time evolution of a single 
peak for a fixed background H2 photodissociating flux and 
also how the internal properties of a given cloud change 
at fixed redshift when the level of H2 photodissociating 
flux is varied. Armed with these examples of the inter- 
nal evolution and structure of the collapsing clouds, we 
discuss in Section §^ the neglect of self-shielding and 
photodetachment in our simulations and then use a simple 
analytical argument to illuminate the essential physics of 
the collapse. We summarize our results in Section §|^. 

2. THE SIMULATIONS 

A self-consistent three dimensional cosmological simu- 
lation of the collapse of low mass (10^ - IO'^Mq) proto- 
galaxies at high redshift (30 > 2 > 19) requires both large 
enough simulation volume to model the gravitational tidal 
forces at work during the formation epoch and high spatial 
(and mass) resolution to probe within the small collaps- 
ing clouds. One technique that has recently been able to 
achieve the large dynamic range necessary to shed light on 
the collapse and fragmentation of these first star produc- 
ing structures is adaptive mesh refinement (AMR). The 
three dimensional AMR algorithm used in our simulations 
is described more fully elsewhere (Bryan 1999; Bryan & 
Norman 1997, 1999; Norman & Bryan 1998). Briefly, the 
code uses an adaptive hierarchy of rectangular grid patches 
at various levels of resolution very similar to the algorithm 
described by Berger & CoUcla (1989). Each grid patch cov- 
ers some region within its parent grid needing additional 
refinement and may itself become a parent grid to an even 
higher resolution child grid. We take the ratio of parent 
to child grid mesh spacing to be two. The dark matter 
is followed using methods similar to those presented by 
Couchman (1991) and the gas hydrodynamics uses the ar- 
tificial viscosity method of Stone & Norman (1992). 

We use a flat, low matter density ACDM model for 
our simulations with ilo = 0.4, fib = 0.05, JIa = 0.6, 
h = 0.65, (Tg = 0.8, and n — 1 where Qq, fib, and flA are 
the fraction of the critical energy density carried in non- 
relativistic matter, baryons, and vacuum energy, respec- 
tively, h is the dimensionless Hubble parameter in units of 
100 kms~^ Mpc""'^ , (78 is the density fluctuation normal- 
ization in a sphere of radius 8h~^ Mpc, and n is the slope 
of the primordial density perturbation power spectrum. 
The parameters are chosen to provide good consistency 
with observation. The fluctuation normalization is con- 
sistent with the CMB quadrupole as measured by COBE 
(Bunn & White 1997) and also with observations of the 
number density of galaxy clusters (White, Efstathiou & 
Frenk 1993; Bond & Myers 1996). f^^/i^ is consistent with 
big bang nucleosynthesis and the measured abundance of 
primordial deuterium (Copi, Schramm, & Turner 1995; 
Buries & Tytler 1998), and JIa is consistent with the up- 
per hmit (17a < 0.7) of Maoz & Rix (1993) and the best 
fit parameters of Ostriker & Steinhardt (1995). 

The simulations are initialized at z = 99 with density 
perturbations generated for the above ACDM model using 
the Eisenstein & Hu (1998) transfer functions in a comov- 
ing simulation volume of 1 Mpc'^ . This simulation volume 
is large enough to ensure that the fundamental mode in 



the box remains linear at the lowest redshifts considered 
here. A low resolution run is used first to identify a re- 
gion of active structure formation. In order to achieve 
high mass resolution in the region of interest the simula- 
tion is reinitialized with multiple static refinement levels 
surrounding the chosen region such that each successive 
static parent grid contains the lagrangian volume of its 
child grid. This results in a mass resolution in the initial 
conditions for the region of interest of 4.78(38.25)M0 for 
the gas (dark matter), respectively. In addition, that re- 
gion is allowed to dynamically refine so that the local Jeans 
length is resolved at all times by at least four grid zones 
and no grid cell contains more than four times the initial 
gas mass element (4.78Mq) or ten times the initial dark 
matter element (382. 5M0). We limit the total refinement 
to 14 levels within a 64^^ top grid which results in a maxi- 
mum dynamic range of 1, 048, 576. The comoving spatial 
resolution of 0.95 pc at maximum refinement translates 
into a physical spatial resolution of 0.03 pc (0.05 pc) at 
z = 30 {z — 19), respectively. We call a peak "collapsed" 
when it reaches maximal refinement. Once this occurs, 
we still wish to follow the subsequent evolution of struc- 
ture in the region, but are uninterested in the details of 
the collapsed peak. In order to prevent the peak from col- 
lapsing further (which could cause numerical instabilities), 
we introduce a form of artificial pressure support. This is 
done by defining, for each cell, an effective pressure which 
is the greater of the thermal pressure and KGplAx'j/fi, 
where pb is the baryon density in the cell. Ax/ is the cell 
width on the finest level and n = 1.22mH is the usual 
mean mass per particle. The dimensionless constant K is 
set to 100, a value which spreads the mass over a spheri- 
cal region with a radius of several cells. The form of this 
expression is chosen by matching the thermal energy and 
gravitational self-energy of the gas; notice also that it re- 
sults in a polytropic equation of state. We also smooth 
the gravitational potential by 4Aa; / which helps to reduce 
the amplitude of the artificial pressure support required. 
Although introducing artificial pressure support causes us 
to lose information about the detailed morphology of the 
inner 1 - 2 pc (~ 1 % of the virial radius) of the peak 
after collapse, it should not affect the determination of 
cooled gas fractions in Section §|| since the cooling region 
lies outside of the collapsed H2 core. 

Our data set consists of results from four simulations 
starting from identical cosmological initial conditions on 
the density fields. Three of the simulations follow the 
nonequilibrium, time dependent evolution of nine chem- 
ical species (H, H+, He, He+, He++, e~, H2 , H+, H") 
using the algorithm of Anninos, et al. (1997) initialized 
with post-recombination abundances (Anninos & Norman 
1996). The reaction network for even this simple chemistry 
is complex. Numerous authors have worked to identify 
the dominant reactions and calculate their rate coefficients 
including recently Haiman, Thoul & Loeb (1996), Abel 
et al. (1997) and Galh & Palla (1998). We use primarily 
the rate coefficients of Abel et al. (1997) with the excep- 
tion of rate coefficients for the processes H + e~ ^ H~ +j, 
i?2 + e" ^ 2H + g-, and H' + H+ + e" which 

are taken from Shapiro & Kang (1987). We have also in- 
cluded the three-body process 2H + H ^ H2 + H (Orel 
1987, Palli et al. 1983), although its contribution is neg- 
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ligible at the densities considered in this work. The small 
post-recombination abundance of H2 in the intergalactic 
medium (~ 2 x 10~^) is too low for cooling to be effec- 
tive. However, at the densities and temperatures in the 
interiors of collapsing proto-galactic clouds, H2 formation 
is catalyzed by ions through the process 

H + e-^H-+-f (1) 

H + H-^H2 + e~. (2) 

and can reach abundances high enough to act as an effi- 
cient coolant allowing the first stars to form. 

Once the first stars form, radiation above 13.6 eV pro- 
duced by these stars is quickly absorbed by the dense 
neutral hydrogen surrounding them. However, the inter- 
galactic medium is optically thin to radiation below the 
HI ionization limit. This radiation escapes to form a soft 
ultraviolet background of radiation throughout neighbor- 
ing regions with sufficient energy to photodissociate H2 by 
means of the Solomon process 

H2+j^H;^H + H, (3) 

thus inhibiting continued star formation (Haiman, Abel & 
Rees 1999). H2 is any of 76 Lyman- Werner resonances in 
the energy range 11.18 ~ 13.6 eV . We use a very simple 
model to investigate the negative feedback effects of this 
soft UV radiation. The radiation field is assumed to be 
constant throughout the simulation with mean photon en- 
ergy of 12.87 eV . In terms of the constant mean flux F^w , 
the H2 photodissociation rate coefficient for the Solomon 
process denoted in Equation || is (Abel, et al. 1997) 

fcrf,,, = 1.1 X 10*^FwS-i. (4) 

We distinguish three of our simulations by different lev- 
els of this background flux, i.e. Flw = (ffist structure 
formation), 10^^^, and 10~^^ ergs"^ cm~^ Hz~^ , respec- 
tively. These flux levels correspond to J21 normalizations 
of the spectral intensity (Haiman, Abel, & Rees 1999) of 
J21 = 0, 10-2-1 and lO-^-^ ergs'^ cm'^ Hz"! sr"! , re- 
spectively. We also include results from a fourth simula- 
tion that starts from the same cosmological initial con- 
ditions, but evolves only six species (H, H+, He, He+, 
He"*""*", e~) representing the extreme case when destruc- 
tion of molecular hydrogen is complete. 

We neglect the time dependent turn on and build-up 
of the soft radiative field. Thus the UV radiation experi- 
enced by the collapsing structures may be interpreted as 
due either to "direct" , nearby sources or a more diffuse soft 
background. We also neglect the effect of self-shielding by 
the collapsing cores. As will be discussed more fully in 
Section §||, this should have only a small affect on our re- 
sults. In any case each of these neglected effects would 
tend to reduce the impact of the photodissociating flux on 
the collapsing clouds. Thus our results for the simulations 
with Flw = 10~22 and 10"^^ erg s~-^ cm"^ Hz~^ should 
be interpreted as representing the maximal suppression of 
structure formation at these flux levels. 

We also neglect the photo-detachment of H~ by the 
UV radiation field. One might argue that the presence 
of the radiative background might photo-detach H~ be- 
fore it could catalyze H2 formation by the process given 



in Equation g. This could hinder the collapse of clouds 
by preventing the build-up of coolant in the protogalac- 
tic cores. While the photons clearly have sufficient energy 
{Ej > 0.755 eV ) to photo-detach H~, we will show in 
Section §^ that H~ photo-detachment is suppressed by at 
least two orders of magnitude relative to H2 formation at 
the temperatures, densities and UV background flux levels 
we consider. 

3. THE SAMPLE OF SIMULATED HALOS 

The locations of high density pregalactic clouds are iden- 
tified in the simulations using the HOP algorithm (Eisen- 
stein & Hut 1998) acting on the dark matter density 
distribution. We find that the high density peaks form 
within and at the intersection of sheets and filaments as 
expected in any cold dark matter dominated cosmology 
(Abel, et al. 1998). We choose the 8 -10 most massive 
peaks of the ~ 200 identified by HOP at each redshift for 
further analysis. In Figure |l| we present the mass of these 
clouds as a function of redshift for the simulation with- 
out H2 . Since we use identical initial conditions for each 
of the four simulations, this is essentially the same set of 
peaks used for analysis in the simulations with H2 present 
and varying levels of the radiative background flux. These 
peaks span a mass range from 3 x 10'* to 7 x IO^Mq and 
are modeled in our simulations by ^ 1000 to 125, 000 dark 
matter particles, respectively. The absence of low mass 
peaks in the bottom portion of Figure ^ is an artifact of 
our selection criterion. However, the increase in the max- 
imum mass of the cloud with decreasing redshift is real 
and characteristic of the "bottom up" growth of structure 
in these models. Small clouds form first at high redshift, 
merging to form more massive structures later. There is 
clear evidence for mergers among the clouds we consider. 
See, for example, the merger of the two most massive peaks 
in Figure |l| into a single peak between redshifts z = 21 and 
z — 20.5. However since we do not follow all of the low 
mass clouds in the simulation, we do not construct com- 
plete merger trees for the sample. Thus we caution the 
reader that all the peaks are not statistically independent 
of each other. This is particularly noticeable for the high 
mass peaks which are rare in our simulation volume and 
may result in a decrease in the scatter of measured prop- 
erties for the higher mass clouds. 

We define the virial radius, ryir, as the radius of a sphere 
within which the average density in the cloud exceeds 200 
times the mean matter density of the universe. The virial 
mass, Mmr is the total mass contained within that ra- 
dius. Before cooling becomes important the mean gas- 
mass- weighted temperature of the cloud obeys a virial re- 
lationship between temperature, mass and redshift simi- 
lar to that observed in galaxy clusters (Bryan & Norman 
1998). For a neutral mix of primordial gas 

T.„(K ) = l55/(/^2A,i?2)V3(_^)2/3 (5) 

where / = 0.8 is a normalization factor, h = 0.65 is the 
dimensionless Hubble parameter at z = 0, Ag = 200 is the 
collapse overdensity and 

E^ ^ no{l + zf + Ha ^ no{l + zf (6) 
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for no = 0.4 and f^A = 0.6 at redshifts 30 > z > 19. 
In Figure ^ we compare the mean gas-mass- weighted tem- 
perature Tmean, for the cfouds in our sample to the virial 
temperature T^j^ defined by Equation ^ for the case when 
no molecular hydrogen is present and for the case F^w — 
when molecular hydrogen cooling is maximal. All of the 
clouds in our sample have virial temperatures Tyir < 8000 
- 10, 000 K , the temperature necessary for HI line cooling 
to become important. Peaks analyzed from the simulation 
without II2 can not cool and their mean temperatures are 
well characterized by the virial relation given in Equation 
1^. Once cooling becomes important, the mean gas- mass- 
weighted temperature of the cloud will drop significantly 
below its virial value as demonstrated in Figure by peaks 
from the simulation that includes the H2 chemistry. 

In Table we compare the properties of the first pro- 
togalaxies to "collapse" , i.e reach maximal refinement, for 
the simulations with II2 cooling and soft UV background 
flux levels Flw = 0, 10^22^ ^^^^j ^g^^i ergs"^ cm'^ Hz^i . 
Protogalactic clouds in the simulation without H2 have not 
collapsed by redshift z = 19. The first structures {Flw = 
0) in this ACDM model collapse early (z ^ 30) and are 
smaU {ryir ~ 50 pc), low mass {Myir ~ 1.4 x IO^Mq) 
clouds. Simulations of the minimum collapse mass by 
Fuller & Couchman (2000) find the minimum mass to col- 
lapse at z = 30 to be 3 X 1O^M0 (expected to be robust 
within about a factor of two) thus consistent with our re- 
sults. Abel, Bryan, & Norman (2000) find a 7 x IO^Mq 
peak collapses first (at z = 19) in their high resolution 
SCDM simulation. This too is consistent with the mini- 
mum mass collapse curves of Fuller & Couchman. Thus de- 
spite differences in numerical method and cosmology there 
appears to be general agreement that the minimum mass 
of the first star forming structures is ^ 1O^M0. This is 
about an order of magnitude below the earlier estimates of 
Tegmark et al. (1997) due primarily to the use of improved 
II2 cooling functions. 

As the amount of the H2 photodissociating flux in- 
creases, the mass threshold for collapse increases pushing 
the redshift of collapse to lower values. However, even with 
flux levels as high as Flw = 10^^^ ergs"^ cm~^ Hz~^ , 
H2 can survive in the cores of the clouds such that ob- 
jects {My^r ^ 1.26 X IO^Mq, T„r ^ 2262 K ) below the 
mass threshold for HI line cooling are still able to col- 
lapse. However, the collapse in this case is delayed until 
z = 21.5. This is consistent with the 1-D simulations of 
Haiman, Abel & Rees (1999). See, in particular, their 
Figure 6 for J21 — 10^^'^. Furthermore the ability of a 
cloud to collapse does not seem to be strongly dependent 
on the angular momentum L of its dark matter halo, since 
the threshold objects we found to collapse spanned the 
fufl range of expected spin parameters 0.02 < A^m < 0.11 
where A = {L\E\^/^)/{GM'^/^), E is the total energy, M 
is the mass, and G is the gravitational constant. 

In order to check the numerical robustness of our re- 
sults, we ran several lower mass resolution simulations 
(AM = 306 (38)M0 for dark matter(gas)) using the same 
random seed for the initial conditions and region for dy- 
namical refinement as the high mass resolution runs pre- 
sented here. We also ran several different realizations 
of the initial conditions at this lower resolution for the 
Flw = 10~^^ ergs""'^ cm"^ Hz^^ case to ensure that there 
was nothing exceptional about the particular realization 



we present. We found that, although the properties of the 
first peak to collapse were very similar for the two mass 
resolutions, the higher mass resolution run delayed the red- 
shift of first collapse to slightly lower values (Az ~ 1). 
This difference in the redshift at which maximal refine- 
ment is reached, however, is already less than the cosmic 
variance (Az ^ 2) in this quantity observed between sim- 
ulations run with different realizations of the initial con- 
ditions. The mass threshold for collapse was also reduced 
by « 30 % in the higher mass resolution run compared 
to the lower resolution simulation. This is slightly larger 
than the ±20 % variation in mass found from simulations 
with fixed mass resolutions but initial conditions generated 
using different random seeds. 

4. COLD GAS FRACTION 

A critical input parameter into semi-analytical models 
that seek to include stellar as well as radiative feedback 
is the amount of gas in these first pregalactic objects that 
is available for star formation. Abel et al. (1998) esti- 
mated the fraction of gas able to cool by adding up all 
the mass within the radius of a 4cr peak (in three different 
simulations) for which the temperature in radially aver- 
aged profiles began to decrease. However, cooling alone is 
not sufficient to determine whether the gas is available to 
form stars. It must also become dense. We improve upon 
these previous estimates by measuring both the fraction of 
gas that is able to cool via molecular hydrogen cooling, /c, 
and the fraction of gas that becomes both cold and dense, 
fed, hi the clouds comprising the data sample described in 
Section It is this latter fraction fed of cold, dense gas 
that is available to form stars in these objects. 

In Figure || we present both these statistics as functions 
of the cloud mass M^r for various levels of the background 
radiative fiux Flw- We define the following: 

• fc is the fraction of total gas within the virial ra- 
dius with temperature T < O.bTyir and gas density 
p > lOOOpmean whcrc Pmean IS the mean gas density 
of the universe. This is the amount of gas within the 
cloud that has been able to cool due to molecular 
hydrogen cooling. 

• fed is the fraction of total gas within the virial ra- 
dius with temperature T < O.bTyir and gas density 
p > W^^Mq Mpc-3 - 330 cm-3 . This is the frac- 
tion of gas within the cloud that is available for star 
formation. 

The above criteria to determine the amount of gas in 
each category (cooled or both cold and dense) are applied 
on a cell by cell basis within the virial radius of each peak 
in the data sample. The temperature threshold has been 
chosen to ensure that the gas is substantially cooler than 
the virial temperature of the peak as given in Equation ^ 
The density threshold for cold, dense gas corresponds to 
the gas density at which the baryons become important to 
the gravitational potential and thus to the subsequent evo- 
lution of the core. However, setting a density threshold for 
gas that has been able to cool through the H2 chemistry 
but may not be dense is more subtle. We must be care- 
ful to exclude gas that is cool only because it belongs to 
lower mass clouds that are infalling onto the more massive 
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structure. With the high mass resolution in our simula- 
tions, this substructure due to merging is resolvable. Since 
there is no cooling in the peaks extracted from the simu- 
lation without H2 , any cold gas present must be due to 
infall and merging in the outer regions of the cloud. We 
confirm that this is indeed true in the data set without 
H2 cooling by looking at radial profiles for the cold gas 
fraction in several peaks. We then used this data to deter- 
mine the density threshold for cooled gas given above that 
first minimizes the cold infalling component. As shown in 
the bottom panel of Figure ||, the resultant background 
due to merging is negligible for all but the most massive 
peaks in the sample. The bump in fc in the high mass 
end of the bottom panel is due to the merger of the two 
most massive peaks in the box at z ~ 20.5 seen in Figure 
1^. Even in this extreme case the above density threshold 
keeps the infalling cold component of fc limited to a few 
percent. Note also that fed, the fraction of gas that can 
become dense and form stars, is zero for all peaks in the 
simulation without H2 , as expected. 

Figure ^ represents one of the key results of this paper. 
Both the fraction of gas fc that can cool and the fraction 
fed of cold, dense gas available for star formation increase 
logarithmically with the pregalactic cloud mass over the 
mass range dominated by H2 cooling. This dependence 
can be summarized in a simple fitting formula: 

/, = B, HM/Mth) (7) 

for i ~ c (cooled gas) 01 i — cd (cold, dense gas) and 
M > Mth- The best fit slopes Be for the fraction of 
gas fc that cools are 0.138 ± 0.006, 0.139 ± 0.008, and 
0.085 ± 0.009 for radiative background fluxes of F^w = 0, 
10^^^, and 10^^^ ergs"^ cm"^ Hz"^ , respectively, where 
the quoted errors are standard errors of the fit. How- 
ever, the fraction of cold, dense gas fed increases much 
less rapidly with mass. Best fit slopes Bed for this cold, 
dense component are 0.058 ± 0.006, 0.066 ± 0.007, and 
0.030 ± 0.004 for radiative background fluxes Flw — 0, 
10"^^, and 10~^^ ergs~^ cm~^ Hz~-^ , respectively. The 
equality of the slopes for flux levels Jlw — and 10~^^ 
erg s~^ cm~^ Hz~^ is striking and probably reflects the fact 
that the cooling chemistry is the same in each case in- 
dependent of radiative flux level. The apparent soften- 
ing of the slopes for F^w = 10~^^ ergs~^ cm~^ Hz~^ is 
most likely an artifact of the small number of independent 
high mass clouds in our simulation volume. The mass 
thresholds Mth for nonzero fc and fed are about equal 
at fixed Flw and increase as the level of soft UV flux is 
increased. This is a characterization of the negative feed- 
back expected due to the photodissociation of H2 in the 
lower mass clouds. Our results for the mass threshold can 
be parameterized as; 

MTHiMQ) = 1.25 X 10^ + 8.7 x W^^^f -^' (8) 

The presence of a soft H2 photodissociating flux at the 
levels considered here delays, but does not prevent gas in 
these pregalactic objects from becoming dense due to H2 
cooling in their cores. The fraction of gas available for star 
formation due to H2 cooling can be represented as: 

fcd-OMHM/MTH) (9) 
with AIth given by Equation 



5. UNDERSTANDING PEAK FORMATION AND COLLAPSE 

In this section we use spherically averaged radial pro- 
files of the pregalactic cloud properties to illustrate the 
internal dynamics of its collapse. We first consider the 
evolution with redshift of a collapsing pregalactic object 
when exposed to a fixed level of photodissociative fiux. We 
then consider how the internal properties of a given cloud 
change at fixed redshift when the level of the photodisso- 
ciating flux changes. 

5.1. Evolution with redshift 

In three of the panels of Figure §0 we show the evolu- 
tion of the gas density, temperature, and H2 mass fraction 
from z = 30 to maximal refinement at z = 21.5 for the first 
object to collapse at a radiative flux level Flw — 10~^^ 
ergs~^ cm" -^Hz-i . (See Table |l|). In the fourth panel 
(lower right) we show the evolution of the cooling and dy- 
namical timescales for 25 < z < 21.5 after cooling begins. 
Starting at large radius the upper set of curves represent 
the H2 cooling time teooi given by 

1.5n„fcRr 

teool = , ' 10 

AnHnH2 

where n^, nn, and nH2 are number densities for the gas, 
neutral hydrogen, and II2 respectively, fc^ is the Boltz- 
mann constant, T is the temperature, and A is the H2 
cooling function. The lower set of curves at large radii 
show the dynamical timescale tdyn where 

for the same redshifts. The solid (dot-dashed) horizontal 
line is the Hubble time 

tHubble = T^iW + zfy'^^ (12) 

for z = 21.5(25), respectively, and the heavy dashed line 
is the difference between the two, i.e. the elapsed time be- 
tween the onset of cooling and collapse. The vertical line 
locates the virial radius ryir = 141 pc at z = 21.5 when the 
structure has grown to a mass of 1.26 x IO^Mq. The in- 
creases in density, temperature and H2 mass fraction seen 
outside the virial radius at ~ 1 kpc are due to averaging 
over other massive structures forming nearby. 

During the initial stages of evolution (from 30 < z < 25) 
H2 abundances are too low to be an effective coolant 
and the cloud grows from a mass of ~ 6.5 x 10'* to 
5.3 X IO^Mq through accretion and merging from within 
its host filament. The radial profiles within the core are 
well described by constant densities and temperatures that 
increase as expected as the mass of the object increases, 
while the outer edges of the profiles steepen with de- 
creasing z. The H2 mass fractions are also roughly con- 
stant in the core (mirroring the density profiles) and grow 
with decreasing z. By z = 25 the core temperature has 
risen to about 2900 K , the core gas density has reached 
~ 3 X 10^'^Mq Mpc~3 « 10 cm"3, the H2 mass fraction 



7 



2.3 X 10 ^ and the core has broadened from about 

-ithe 



10 pc to 20 pc. For F^w = 10"^^ ergs^^ cm^^ jjz 
photodissociation timescale 



idiss (-5) 



9 X IQ-y Flw, 



(13) 



is 9 X 10^^ s, much shorter than other timescales in the 
problem. Thus we can understand the H2 mass fraction in 
the core by a simple equilibrium argument. If we assume 
that H2 formation is dominated by the H~ process given 
in Equations |^ and the H2 equilibrium number density 
is given by 



equil 
^H2 



diss 



(14) 



where fcy « 1.8 x io~^^T"-®* cm'^ s~^ over this temperature 
range, and n^r(rie) is the number density of neutral hydro- 
gen (electrons), respectively. At 2; = 25 Equation O pre- 
dicts an equilibrium H2 mass fraction XhT^ ^ -^-^ ^ 10"^ 
in the core that is nearly identical to the value above found 
in our simulation. 

By z — 2A the H2 mass fraction has reached 4 x 10^^ 
and gas in the central region has begun to cool. As cooling 
continues the core density steepens, evolving into a ~ 
profile aX z — 21.5 in the inner 10 pc. The gas density 
reaches number densities n > 4 x 10'^ exceeding the dark 
matter density for r < 1 pc. Notice, however, that most 
of the cooling occurs in the region 2 < r < 20 pc, out- 
side this very dense central core. The temperature profile 
at collapse is very similar to that found by Abel, Bryan 
& Norman (2000) showing with decreasing radius cosmic 
infall, shocking, and a cooling flow. We also find for this 
object that the cooling time becomes comparable to the 
dynamical time over most of the cooling flow region just 
prior to collapse, and that the outer radius of the cool- 
ing region corresponds closely to the radius at which the 
cooling time has dropped below the elapsed time (heavy 
dashed line in Figure 0) in agreement with the collapse 
criterion of Haiman, Abel & Rees (1999). 

Although spherically averaged radial proflles are easy to 
picture, it is important to remind ourselves that the evo- 
lution of these objects is fully three-dimensional and they 
do not collapse in isolation from their surroundings. In 
Figure |^ we show images of the log of the projected dark 
matter density, the baryonic density, and the temperature 
at z — 21.5 for a portion of the dynamically refined region 
in the F^w — 10"^^ erg s~^ cm~^ Hz~^ simulation. This 
region contains the two most massive structures in the sim- 
ulation. The image is 680 pc (proper) on a side and is pro- 
jected along the y-axis. The masses (virial radii) for these 
protogalaxies are 2.4 x 1O^M0(174 pc) for the upper cloud 
and 2.2 x 10*'Mq(170 pc) for the lower cloud, while their 
proper center-to-center separation is only ~ 300 pc. Thus 
their outer regions are already beginning to merge. The 
temperature map shows an intricate pattern of asymmet- 
rical shocks with maximum temperatures around 5400 K 
surrounding elongated knots of cooled gas at the clouds' 
centers. These cooling knots have temperatures of a few 
hundred Kelvin and are the regions of highest H2 con- 
centration and thus gas density. The dark matter ha- 
los, containing ~ 42, 000(37, 000) particles for the upper 
(lower) structures, respectively, clearly show substructure 
near their centers due to the prior infall and merger of 



lower mass clouds that have not yet completely relaxed. 
This merger activity is responsible in part for the aspher- 
ical geometry of the cooling knots and corresponding gas 
density. 

5.2. Changing the photodissociating flux 

In Figure ^ we compare the spherically averaged radial 
profiles of dark matter density, gas density, temperature 
and H2 mass fraction at fixed redshift {z = 21.5) but 
with varying levels of H2 photodissociating flux for the 
2.2 X lO^M© (lower) structure of Figure ^. The cloud 
has only been able to maximally refine in the simula- 
tion with Flw = 0, although it is close to collapse in 
the Flw = 10~^^ case also. This is in spite of the fact 
that the mass of this object is above the masses of the 
first objects to collapse in both the F^w = 10~^^ and 
10~^^ ergs~^ cm~^ Hz~^ simulations (See Table |l|). Prob- 
ably this is due to the strong effect of tidal interactions 
and high level of merging activity in and surrounding this 
object. The features located at 300 pc and 1 kpc in the 
profiles reflect the presence of the upper, 2.4 x IO^Mq pro- 
togalaxy of Figure |^ and the 1.26 x IO^Mq structure we 
analyzed earlier in Figure ^ respectively. The dark mat- 
ter density profiles are essentially the same for all of the 
simulations, independent of H2 fraction or radiative flux. 
This is to be expected as long as the gravitational poten- 
tial is dominated by the dark matter. The dark matter 
density only steepens away from this common proflle in 
the inner '--^ 2 pc of the F^w = case when the gas den- 
sity equals and then exceeds the dark matter density. For 
the Flw = case the gas density steepens in the inner 
several parsecs with gas number densities reaching ~ 10**, 
characteristic of the formation of a dense collapsed core. 
The cooling region extends inward from about 40 pc un- 
til the temperature drops to the minimum value possible 
under H2 cooling, ^ 200 K , at around 3 pc. The H2 mass 
fraction has saturated at ^ 10"'^. The sharp drops in 
temperature that correspond to enhancements in the H2 
mass fraction at 9 and 25 pc are the result of averaging 
over infalling clouds that were able to form H2 in their 
cores and cool before merging with this more massive ob- 
ject. As the radiative flux level is increased, the negative 
feedback becomes apparent. The maximum gas density 
drops in the central region and forms a homogeneous core 
whose core radius also increases as the flux is increased. 
The radius of the cooling region moves inward and cooling 
is suppressed due to the decreased amount of H2 coolant. 
The cloud is clearly in photodissociation equilibrium for 
the high flux (Fr,w — 10~^^) case with the H2 fraction 



given by Equation 14, For the lower flux case, the H2 
fraction is somewhat smaller than the photodissociation 
equilibrium value in the inner 10 pc due to the influence 
of recombination on the formation of H~ ions and thus the 
H2 formation rate. We also see, as expected, that there is 
no cooling at all for the case in which H2 is absent. 

We can gain a better understanding of the above radial 
profiles and their scaling properties by comparing the rel- 
evant physical timescales in the problem. These times are 
shown for both the Flw — 10^^^ case (upper panel) and 
the Flw = 10~^^ case (lower panel) in Figure |7[ In addi- 
tion to profiles for the Hubble, dynamical, photodissocia- 
tion, and cooling times defined earlier, we include profiles 
for the sound crossing time tcross , the recombination time 
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tree given by 



7.7 X 10''r"-''7n, 



and the H2 formation time given by 



tH2 = 



nH2 
kjnHUe 



(15) 



(16) 



All times are in seconds. We also show that the two- 
body relaxation time treiax (Equation 8-71 in Binney & 
Tremaine 1987) for the dark matter particles is long com- 
pared to all other timescales confirming that two-body nu- 
merical effects are unimportant. 

For both cases the photodissociation and H2 formation 
times are much shorter than the Hubble time everywhere 
within the virial radius. In the high flux case they are 
also much shorter than the recombination time so that 
photodissociation equilibrium {tdiss = ^^2) is established 
throughout the cloud (see equation p^). Even so, enough 
H2 forms for the cooling time to drop below the Hub- 
ble time in the inner 40 pc — core collapse is not far 
off. The low flux case is similar in the outer parts of the 
cloud: the recombination time is long and the H2 frac- 
tion is governed by its equilibrium value. However, within 
10 pc recombination becomes faster than photodissocia- 
tion. This reduces the amount of H~ allowed to form 
and suppresses H2 formation below the photodissociation 
equilibrium value. However, enough H2 has formed that 
cooling and collapse proceed faster than photodissociation 
[tdyn < tdiss and teool < tdiss) at these small radii. We 
also find here — as in Abel, Bryan, & Norman (2000) — 
that tcross < tdyn in tbe inner regions of the cloud so that 
the collapse proceeds hydrostatically. 

6. DISCUSSION 

Now that we have detailed examples of the collapse 
properties of clouds in our simulations, we are better able 
to assess the importance of neglected processes in our sim- 
ulations. In this section we discuss in more detail why the 
neglect of self-shielding and H^ photodetachment should 
not significantly change our results. We then use a simple 
analytical model, similar to that by Tegmark et al. (1997), 
to illustrate the essential physics of the collapse. 

6.1. Neglected processes 

We first return to the issue of H2 self-shielding. Al- 
though this effect could, in principle, be important for our 
objects, we argue here that its actual impact is slight. The 
maximum column densities of H2 we encounter in our sim- 
ulations are around lO^'' cm~^, while a column density 
of 5 X 10"'^'' cm~^ is sufficient to start blocking the soft- 
UV flux. So, at first glance, it appears that self-shielding 
should be important; however, this value has been derived 
for a static gas distribution. The Lyman- Werner bands 
consist of a large number of individual lines whose width 
is dominated (in our case) by Doppler broadening. 

There are substantial velocity gradients within all of our 
simulated objects, even those which appear to be relaxed. 
These take the form not just of radial inflow, but also large- 
scale disordered velocity fields (i.e. turbulence). Typical 
rms line-of-sight velocities are around 3-4 km/s. These 
values are larger than the ~2 km/s line widths due to 



thermal Doppler broadening. This means that even if a 
given parcel of gas is behind a large column of molecular 
hydrogen, it may still be photo-dissociated if it has a ve- 
locity shift relative to the blanketing material. There is 
very little continuum opacity. 

It is a difficult radiative transfer calculation to deter- 
mine the precise effect of these velocity shifts; however we 
may consider two limiting cases. The first is the previ- 
ously mentioned static case in which velocities are negli- 
gible compared to the thermal Doppler widths {v/b <^ 1) 
and an H2 column density of 5 x 10^^ cm~^ is sufficient 
to start affecting the photo-dissociation rate. In the other 
extreme, in which v/h 3> 1, the rate is not significantly 
reduced until a substantial fraction of the 11-13 eV range 
is blocked by the Lyman- Werner lines. Since this does 
not occur until the damping wings of the lines become 
important, it requires a much larger molecular hydrogen 
column density, of order 10^^ cm~^ (Draine & Bertoldi 
1996; Glover & Brand 2000). In the simulations, v/b > 1 
so we are likely to be closer to the second regime than the 
first. Even if it were marginally important, there are three 
further reasons to think that self-shielding is unimportant. 
First, these large column densities occur at late times, once 
the core has already collapsed. For most of the evolution 
the objects are clearly in the optically thin regime. Sec- 
ond, the maximum column density quoted above is to the 
center of the object, but the cooling front lies outside of 
the H2 core, and so beyond most of the covering mate- 
rial. Finally, once a sufficient amount of H2 has formed to 
produce a substantial blanketing column, the cooling time 
has dropped below the photo-dissociation time, so that 
photo-dissociation is no longer important (and thus even 
if the soft-UV flux is attenuated, the subsequent evolu- 
tion of the gas would be unchanged) . So while the results 
described earlier are properly lower limits to the cooled 
fraction because of our neglect of self-shielding, a full ra- 
diative transfer calculation should produce nearly identical 
results. 

Let us also reexamine whether H^ photo-detachment 
(which we neglect) could significantly affect the results. 
In order to bracket the size of this effect we consider two 
spectral shapes for the radiation field, J ^ v"^ with a. = 1 
and —1. The softer spectrum, a = 1, is the spectral slope 
expected in this energy range from a typical OB star. The 
steeper slope, a = —1, is the type of spectral shape one 
might expect for a background radiation field built up from 
a number of sources. The cross section (in cm^) for H~ 
photo-detachment is given by (Kang & Shapiro 1987, Abel 
et al. 1997), 



(Jdetaeh = 7.928 X IQ'' [u - uthf'^X^ iov V > 



vth (17) 



where hp is Planck's constant and hpvth — 0.755 eV. The 
H^ photo-detachment rate coefficient kdetaeh can then be 
calculated 



'^detach 



hpv 



<ydetach{v) dv 



(18) 



with hpVH = 13.6 eV. For the two spectral shapes given 
above we find kdetaeh = 1-9 x IQ^^'^F^w for a = 1 
and 6.4 x 1Q~^^Flw for a = — 1 where Flw is the 
radiation flux in ergs~^ cm~^ Hz^^ . 
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The dominant processes other than photodetachnient 
that determine the abundance of H~ are shown in Equa- 
tions |l] and H with rate coefficients kj and k%, respec- 
tively. Furthermore the processes proceed sufficiently 
rapidly compared to the H2 chemistry that uh- is well 
approximated by its equilibrium value (Abel et al. 1997) 
such that (including photo-detachment): 



krUHUe 



kfiUH + k, 



detach 



(19) 



We see that photo-detachment can be neglected as 
long as kdetach ^ kgUH- This leads to a condition on the 
hydrogen number density uh ^ kdetach /ks that can be 
checked against the conditions found within the collaps- 
ing structures. For gas temperatures found in our sim- 
ulations, ks w 1.43 X 10"^ cm^s-i (Abel ct al. 1997). 
Thus photo-detachment is unimportant provided that 
uh > 0.045(Fny/10^^^) for the steeper a — —I spec- 
trum or, a much weaker constraint, that nn ^ 1.4 x 
1Q~^{Flw for the softer a = 1 spectrum. Let us 
check these constraints against the hydrogen number den- 
sities found within the first peak to collapse at our high- 
est flux level (Flw = 10~^^ ergs~^ cm~^ Hz~^ ) shown in 
Figure ^. At z = 25, before cooling commences, hydrogen 
number densities in the core region where H2 will form 
already exceed hh ~ 7. This number density satisfies the 
constraint for the steepest spectrum (a — —1) by more 
than two orders of magnitude so that the neglect of H~ 
photo-detachment does not significantly affect our results. 

6.2. An Analytic Estimate of the Collapse Threshold 

Analytical arguments often provide qualitative insight 
into the key physical principles at work in a complex pro- 
cess. Since we have seen that the H2 fraction is close to 
photodissociation equilibrium throughout our objects, we 
use these equilibrium abundances to develop analytical es- 
timates of the critical H2 fraction and mass thresholds for 
collapse. We adopt the simplified collapse criterion, simi- 
lar to that of Tegmark et al. (1997), that an object will 
collapse if the cooling time tcooi , given by Equation hO 



less than the Hubble time t Hubble , 
Thus 



defined in Equation 12| 



LbugkBT 



< 



AnHnH2 3iJ( 



(20) 



The equality gives the critical number density of H2 nec- 
essary for collapse. Solving for n^j* '^^ obtain 



crit _ r, r)r '"'Q^bT 

'■H2 ~ ^-^"J 



UhA 



3/2 



(21) 



In this equation, A is, in our case, the Lepp & ShuU (1983) 
cooling function. In the low density limit {nH2 ^ 10"* 
cm^"^ ) for temperatures 300 K < T < 5000 K of interest 
here it is well approximated by A w 6.0x lO^^^Tfcse^^'^"/"^ 



such that 



„crit 
^H2 



1.3xlO-'he'^'>/^n'J\l + z) 



3/2 



(22) 



If we now equate this critical abundance of II2 required to 
satisfy our collapse criterion to the photodissociation equi- 
librium abundance from Equation |4| and solve for Fj^]^, 



we obtain an expression for the critical flux that just allows 
a cloud with virial temperature T to collapse. 



Tpcrit 



10 



^21 



126n|,a:,r0-88e-730/T 



/ll7y^(l + z)3/2 



(23) 



where the electron fraction — n^/rifj. 

To complete the above relationship we need to approx- 
imate the electron fraction and the gas number density in 
the cloud core. If the collapse occurs well within one re- 
combination time, the electron fraction, Xe, is expected to 
be close to its residual value Xe ~ 1.2 x 10^^r2Q^^/(Ob/i). 
The average baryon number density in the universe is given 
by ns ~ 9.3 x 10"'^(1 -I- zYVlbh"^ cxnT^ . The gas temper- 
ature after thermal decoupling from the CMB at z ~ 200 
but before reionization z > 7 drops adiabatically with 
TiGM ^ 0.0135(1 -I- z)^ K . If we assume that the baryons 
in the central regions of dark matter halos were never 
shocked, we flnd the maximum density that can be reached 
by adiabatic evolution is rimax = nBiTmr/TiGmY^'''^^'^^ ■ 
Here Tyir denotes the virial temperature of the dark matter 
halo and 7 = 5/3 is the ratio of speciflc heats for an ideal 
gas. This gives a collapse redshift independent maximal 
density as a function of virial temperature for adiabatic 
evolution of the cloud as n^ax ^ 1^7 ^Ibh"^ [T /IQQQf/"^ 
cm^^ . We expect adiabatic evolution to be a good approx- 
imation for low mass clouds, but caution that for larger 
mass structures merging and shock heating become im- 
portant. Inserting the expressions for n„iax and Xe into 
Equation 23 we obtain: 



j^crit 



10 



-21 



223nhh 



'(—) 

Vioooy 



3.88 / , ^ ^ -3/2 

g-730/T/_|_ 



20 



(24) 

The structures that are first able to collapse in our sim- 
ulations all have virial temperatures between 700 K and 
3000 K where the exponential in Equation e3 can be ap- 
proximated by a power law, g-'^^o/T ^ 0.48(r/1000)° '^. 
Inserting this approximation into Equation p3 and invert- 
ing, we obtain an analytic approximation for the critical 
virial temperature for collapse at a given soft-UV flux level. 



rpcrit \ 
vir \ 

1000 / 



0.36 {VLbh^ 



Flw 



1 



20 



3/2N 



(25) 

For Flw = 10 ergs ^ cm ^ Hz ^ at z ~ 25 in our 
cosmological model [h — 0.65, Q.b — 0.05) this gives a 
Tyj"* ^ 920 K , somewhat lower than but comparable to 
the Tyir ^ 1470 K found at that redshift for the first peak 
to reach maximal refinement at that flux level in our sim- 
ulation (See Figure ^). The disagreement in detail is not 
surprising since the approximation invokes a simplified col- 
lapse criterion and ignores the nonadiabatic evolution of 
the collapse. However, the analytical argument does serve 
to ellucidate the essential physics at work in the collapse. 

7. SUMMARY 

In this paper, we have used high-resolution numerical 
simulations to investigate the effect of radiative feedback 
on the formation of 10^ — l(f Mq clouds. This range of ob- 
jects is important because they are large enough to form 
molecular hydrogen, but too small to cool by hydrogen 
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line cooling. Therefore they are susceptible to the de- 
struction of molecular hydrogen from soft-UV flux in the 
11-13 eV range (the Lyman-Werner bands). This radia- 
tion comes from the very first generation of stars and so 
its amplitude is unclear; we test four cases, ranging from 
zero flux to complete H2 destruction. The initial condi- 
tions are drawn from a ACDM cosmological model, and 
the simulations evolve the non-equilibrium rate equations 
for twelve species of hydrogen and helium. Our main con- 
clusions may be summarized as follows. 

• A soft-UV flux does delay cooling and star formation 
in the cores of clouds in this mass range. An increas- 
ing amount of radiation delays collapse until later 
redshifts when larger objects have collapsed. How- 
ever, even without any radiative feedback, objects 
with masses less than ~ 1.4 x 10^ Mq (at z ~ 30) do 
not cool efficiently due to their low H2 fraction. We 
provide a fit to the mass threshold for collapse, given 
the mean flux in the Lyman-Werner bands (equation 

• The amount of cold, dense gas which is available for 
star formation in the cores of the collapsed objects 
depends primarily on two numbers, with relatively 
little scatter: the flux of soft UV radiation and the 
mass of the cloud. We show the cold, dense fraction 
increases logarithmically with the mass of the halo 
and we provide a fit for this relation as well (equation 

I). 

• The evolution of these objects is relatively straight- 
forward. With a negligible soft-UV radiation field, 
the formation of H2 continues until it is suppressed 
by the removal of electrons from the core due to re- 
combination. However, by the time this occurs, gen- 
erally enough H2 has formed to allow core collapse. 
When radiative feedback is important {Flw ^ 5 x 
lO^^'' ergs~^ cm~^ Hz^^ ), the destruction and for- 
mation processes are sufficiently fast so that the H2 



fraction is in equilibrium (equation |IJ) and recom- 
bination plays little role until after core collapse be- 
gins. This allows the derivation of an analytic esti- 
mate for the collapse threshold which comes close to 
matching the full simulation results. 

Finally, there are a number of other results which de- 
serve further investigation. One of these is that the col- 
lapse redshift does not seem to depend on the angular mo- 
mentum of the halo (see section Another is that we 
do not appear to witness fragmentation in these objects. 
The few times in which we clearly see more than one col- 
lapsed core in a halo appears to be due to the merging of 
a number of halos, each of which already contained a core. 

We also examined a number of process which we omit- 
ted from the simulation, including H2 self-shielding and 
photo-detachment, none of which appears to be impor- 
tant for this calculation. These results were derived within 
the context of a currently popular A-dominated cosmologi- 
cal model. However, at these high redshifts, the effective $7 
should be close to one, so the most important parameter is 
the amplitude of the power spectrum on these scales. Rela- 
tively small changes to the amplitude (such as would arise 
from the current uncertainty in erg) are unlikely to have 
a major impact on the results presented here, although 
it should always be kept in mind that we are extrapolat- 
ing the power spectrum to scales significantly smaller than 
constrained by observations. 
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Table 1 

Properties of the first object to reach maximal refinement for various levels of soft-UV 

BACKGROUND FLUX Flw- 



-F'LwCerg/s/cm^/Hz) 


2 


Mtotal 












fed 







30.2 


1.37 X 10^ 


713 


49 


1.28 X 10^ 


8800 


616 


0.07 


0.023 


10-22 


24.4 


8.83 X 10^ 


2016 


111 


8.02 X 10^ 


8.08 X 10"' 


1734 


0.021 


0.107 


10-21 


21.5 


1.26 X lO'^ 


2262 


141 


1.15 X lO'* 


1.15 X 10^ 


4180 


0.037 


0.042 



Note. — z is the redshift of maximal refinement. T^ir is tlie virial temperature given by Equation pi Mtotah ^dmy and Mgas 
are the total, dark matter, and gas masses, respectively, in Mq within the virial radius r^ir-. M^d (jWq) is the total amount of 
cold, dense gas (T < 0.5T^ir, p > 10^^ Mq Mpc— ^) within the virial radius, fed is the corresponding cold, dense gas fraction, 
and Adm is the dark matter halo spin parameter. 
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Fig. 1. — The virial mass as a function of redshift for peaks used in the analysis sample for the case with no H2 present. 
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Fig. 3. — Fraction of cold gas within the virial radius as functions of cloud mass and soft-UV background flux F^^y in ergs~^ cm~^ Hz"-' . 
Open symbols represent /c, the fraction of gas that has cooled via H2 cooling (T < O.bTyir, p > lOOOpmean with pmean the mean baryonic 
density of the universe). Filled symbols represent fed, the fraction of cold, dense gas (T < O.bT^ir, p > W^^ Mq Mpc~^) available for star 
formation. 
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Fig. 4. — Evolution of the radial profiles for the first peak to collapse at radiative flux level Fz,vk = 10 ergs ^ cm" -2 Hz" 1 from z = 30 
to maximal refinement at z = 21.5. The heavy horizontal lines in the lower right panel represent the Hubble time at z = 25 (dot-dash) and 
at 2 = 21.5 (solid), while the heavy long-dashed line is the elapsed time from the onset of cooling until collapse. The vertical line marks the 
virial radius at z = 21.5. 
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Fig. 6. — Radial profiles of the daj^ matter density, gas density, temperature and H2 mass fraction for the lower (M^ir = 2.2 X IO'^Mq, 
^vir = 170 pc) structure of Figure B for various levels of the H2 photodissociating flux F]^\y in orgs^^ cm"^ Hz~^ . The photodissocia- 
tion equilibrium values for the H2 mass fractions are also shown for Fl\y = 10"^'^ ergs~^ cm^^ Hz^^ (open circles) and Fj^w = 10^^^ 
ergs~^ cm^^ Hz~^ (filled circles) for comparison. 
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